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Abstract 

The morphology of adiabatic accretion flows with small non-zero specific angular 
momentum has been investigated in the axisymmetric and nonviscous limit. For an 
initial state characterized by a Bondi flow with the specific angular momentum dis- 
tributed with respect to polar angle, a travelling shock wave forms which propagates 
more rapidly in the equatorial plane than in the plane perpendicular to it resulting 
in the formation of a hot torus. In cases where the incoming flow is restricted to lie 
near the equatorial plane, a strong wind forms directed away from this plane with 
the tendency for the formation of a nonsteady shock structure. As the height of the 
incoming flow is increased it is found that the resultant wind weakens. The parame- 
ter regime which delineates accretion flows characterized by a travelling shock and a 
nearly standing shock structure is presented. 

keywords: accretion, accretion disks — hydrodynamics — methods: numerical 
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1 Introduction 

The process of accretion is central to our understanding of many diverse astronom- 
ical objects and phenomena. The formation of young stellar objects on the galactic 
scale and the underlying cause of the activity of galactic nuclei on the extragalactic 
scale are just a few such examples. A realistic description of accretion flows is likely to 
be multi-dimensional and time dependent in character and may involve the interplay 
between various complex physical processes such as dissipative magnetohydrodynam- 
ics and radiation. However, in some idealized cases simple geometry or physics may 
suffice to facilitate solution by analytical or semi-analytical means. Among such prob- 
lems include spherical accretion in the adiabatic or isothermal approximation, and 
disk accretion driven by viscous effects. 

The classical case of spherical accretion onto a stationary point mass in the hy- 
drodynamical approximation was first considered by Bondi (1952). In this study he 
showed the existence of several families of solutions and described their topology. Of 
import is the transonic solution which bears his name. The extension of the accretion 
problem to the two dimensional case in which the accretor is in relative motion with 
respect to the surrounding medium was actually investigated earlier in the particle 
approach by Hoyle & Lyttleton (1939) and Bondi & Hoyle (1944). Later, Ruderman 
& Spiegel (1971), Wolfson (1977) and Bisnovatyi-Kogan et al. (1979) reexamined 
this problem semi-analytically in the hydrodynamical approximation. However, no 
analytical solution exists in the general case in which the hydrodynamical effects are 
fully considered. In fact, the recent numerical calculations reveal solutions which are 
complex and time dependent. In particular, the flows are found to exhibit a wide 
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range of structures including shocks, jets, and transient accretion disks (see, e.g., 
Shima et al. 1985; Taam & Pryxell 1988; Ruffert & Arnett 1994). 

Accretion flows with large speciflc angular momentum have also been the subject 
of intensive theoretical study following the seminal analytical work of Shakura & 
Sunyaev (1973). In the case where matter is viscously driven recent numerical studies 
have focused on the formation of a disk wind (Eggum, Coroniti, & Katz 1987, 1988) 
and meridional flows in geometrically thin accretion disks (see Kley & Lin 1992; 
Rozyczka, Bodenheimer, & Bell 1994). In such studies viscosity and radiation transfer 
are essential for a proper description of the disk. On the other hand, if the speciflc 
angular momentum of the accreting matter is low, then the flow is, in general, pressure 
driven and the disk is geometrically thick. Accretion in this regime was studied in the 
pioneering investigations by Hawley, Smarr, & Wilson (1984a,b) and Clarke, Karpik 
& Henriksen (1985), and recently by Molteni, Lanzafame, & Chakrabarti (1994), 
Ryu et al (1995), and Igumenshchev, Chen, & Abramowicz (1996). The appearance 
of vortices and shock waves are common features in these accretion flows. 

The results from such investigations, taken as an aggregate, do not present a clear 
picture of disk accretion in the geometrically thick disk limit. For example, Hawley 
et al. (1984b) studied the time dependent two dimensional flows in the vicinity of a 
black hole in the limit that the speciflc angular momentum of the accreting matter lies 
between that of the marginally stable orbit and the marginally bound orbit. In this 
case, a shock front forms and travels outward in the disk. This result was conflrmed in 
an independent study by Clarke et al. (1985) in the Newtonian approximation. Their 
results, however, are in direct contrast to the results obtained in the recent studies 
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of Molteni et al. (1994) and Ryu et al. (1995). In both these latter investigations, 
a standing shock front formed in the disk confirming the earher work of Chakrabarti 
(1990) and Chakrabarti & Molteni (1993). Furthermore, both studies demonstrated 
that a wind emanating from the inner disk region is formed. Further complicating 
the description of the accretion process is the recent study by Igumenshchev et al. 
(1996) who investigated similar thick disks, but with viscous effects included. In this 
case shock formation was absent in the accretion flow in conflict with the conclusion 
reached in Molteni et al. (1994). 

In this paper we report on the results of two dimensional hydrodynamical simu- 
lations which delineate the parameter regime where accretion flows are characterized 
by a travelling shock front and a standing shock front in the inner disk. To make 
comparisons more meaningful to previous studies, we restrict our study to accretion 
flows with small non-zero speciflc angular momentum in the adiabatic, nonviscous 
hydrodynamical approximation. Our goal is to provide for an understanding of the 
various accretion flows in terms of the vertical extent, specific angular momentum 
of the incoming matter (see Hawley et al. 1984a,b), and the boundary conditions 
imposed on the flow. In §2 we outline the speciflc approximations underlying our 
investigation. The initial conditions and the boundary conditions for the suite of 
problems is outlined in §3 and their numerical results are presented in §4. Finally, we 
discuss our results in the context of previous studies in the last section. 

2 Formulation 
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In this investigation it is most convenient to use spherical coordinates {r,9,(f). 
The accretion disk is assumed to be non self-gravitating and axisymmetric with the 
rotation axis coincident with the polar axis (6* = 0). The flow is assumed to be 
adiabatic with no local radiative coohng included. The gravitational field of the 
black hole is described in terms of a Newtonian potential in order to compare to 
previous studies 



where M the mass of the central object. 

The governing equations describing the flow in these coordinates can be expressed 
in the following form. For the equation of continuity, 

dp I d , r, ^ Id, , , 

-£ + ^^^V + —^^{pvesme) = 0. 2.2 
at or rsmOoO 

The equations of motion take the form 

d , s I d , n Id, . dp pive^ + vJ) ,^ 

-^ip^r) + -^^{r^pvD + — — {fyvrvesme) = -p— - / + ^, 2.3 

at r^ or r sm 9 oO Or Or r 

^ypve) + pvevr) + —^^{pv^sme) = — — \ ^ , 2.4) 

at r^ or r sm 9 o9 r 09 r r 

i(P^) + \^i^'P^^r) + -^^{pivesm9) = 0. (2.5) 
at r^ or r sm 9 09 

The energy equation is written as 

d Id Id 

^(P^) + — ^(^Vt^r) + ^^^(p£^;esin^) = -pdivy, (2.6) 
at r^ or r sm 9 o9 
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with the equation of state taken to be, 

p=(7-l)p£. (2.7) 

Here p, {vr, vq, v^), e, and p are the density, the velocities, the specific internal energy, 
and the pressure of the flow respectively. In addition, i = Vi^r sin 9 is the specific 
angular momentum, and 7 is the adiabatic index assumed to be a constant (7 = 1.5 
or 5/3). 

For definiteness we have fixed M — IOMq and i{rout: 7''/2) — io — '2.Rgc and have 
used Rg and Rq/c for the units of length and time respectively. Here, Rq = 2GM/c^, 
is the Schwarzschild radius for a black hole of mass M. In these units, the velocity 
is expressed in terms of the speed of light, c, and the energy per unit mass in terms 
of c^. For variations in M and £0, the results can be scaled by changing the units of 
length and time to (^o/2-Rgc)^-Rg and {Iq/2RqcYRq/ c respectively. 

The above equations are solved using the explicit second-order Eulerian hydro- 
dynamical PPM method (Woodward & Colella 1984). In the present study, the 
calculational domain is located within rj„ < r < Tout a-nd < 9 < Tr/2, where rj„ and 
rout define the inner and outer boundaries respectively. The numerical grid in the 9 
direction is homogeneous, while in the radial direction, the width of the radial zone 
is given by Ar = rA9/{l — 0.5A9). We assume axial symmetry with respect to the 
polar axis ^ = and mirror symmetry with respect to the equatorial plane 9 — 7r/2. 
At the inner radial boundary, zero derivatives with respect to radius r are used for 
variables such as vg and v^, and linear extrapolations are used for p, Vr, and e. At 
the outer boundary, the values of all the variables are fixed in time. 
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3 Initial Conditions 

For a polytropic equation of state, p = Kp"' where K is the polytropic constant, 
the BernouUi function of the flow becomes 

^ = -^Kp(^-'^ + I + I + + (3.1) 

7 — 1 2 2 2(rsin6')^ 

The initial state of the flow can be constructed under additional assumptions. 



3.1 Bondi type conditions 

To construct a Bondi type condition, we first assume ve — and £ — 0. Then 
equations (3.1) and (2.2) become respectively, 

^ -Kp'^-^ + '^-^^B = const, (3.2) 



7-1 2 r 

and 

Anr^pvr = -M, (3.3) 

where M is the mass accretion rate. It is common to assume the sound speed at 
infinity as Coo, which is related to the Bernoulli parameter as = 0^/(7 — 1), where 
the local adiabatic sound speed is = 7Poo/Poo — Kjp^~^\ With one can 
define the Bondi radius as Rb — GM/c^. thus, B/c^ = 2(-^-i) • transonic 
solution which we choose, specification of only two of the three constants, B, M, and 
K is required because the regularity condition at the sonic point provides a constraint 
among the constants which gives the maximal accretion rate. For 7 = 5/3 the sonic 
radius is at the origin. To consider a more general case, we use 7 = 1.5. 



In the case where the specific angular momentum does not vanish, we assume 
that the initial condition can be described by a superposition of the specific angular 
momentum on the Bondi solution. Assuming that vo = and the Mach number, 
the polytropic constant, and the Bernoulli function of the fiow remain unchanged, we 
take the radial velocity and the density to be determined according to equation (3.2). 
Two cases for the distribution of specific angular momentum are considered. In the 
first case 



It should be noted that, in the latter case, when 9 is small, no solution may exist 
since i does not vanish. Under such a situation, Vr is set to zero and a numerically 
small quantity is assigned to p. 

We have also considered cases in which the fiow is restricted near the equatorial 
plane. In these cases, a vertical extent of the fiow, H, is required. For r cos 9 > H, 
the flow is truncated, i.e., the velocities are set to zero and a numerically small value 
is assigned to p. 

3.2 Parallel type conditions 

To compare with the previous study of Ryu et al. (1995), we consider parallel 
type initial conditions in which the velocity of the incoming fiow is parallel to the 




(3.4) 



and the other is described by 



^(r, 9) — Iq — const. 



(3.5) 
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equatorial plane, 



Vr — —V sin 6, vg — —V cos 6, 



(3.6) 



where v — y(v^ + vj) . The Mach number of the flow, M. — y{v^ + vj) / Cg is assumed 
to be constant and is a free parameter to be specified for a solution. Thus, for a given 
polytropic constant, K, solutions of the flow can be solved from a single equation 
(3.1). Here, we have assumed that B — and i — io — const. The vertical height of 
the flow is H. In addition, we have used 7 = 5/3. 

3.3 Torus type conditions 

In the case of a static torus, Vr — and vq — 0. For a given specific angular 
momentum distribution, i = const, equation (3.1) can be solved analytically for 
the density distribution for a fixed B which is a constant here. Specification of the 
solution is then determined by the polytropic constant K. A general property of such 
an initial condition is that solutions cannot exist inside the funnel wall which, under 
the Newtonian potential, is determined by the condition 



In our case, the location of the funnel wall on the equatorial plane is simply ARq for 
B — 0. It should be noted that, in the Newtonian approximation, the funnel wall 

is always closed at the inner edge of the wall independent of B and i 0. Thus, 
accretion is forbidden for a static torus. 



GM 



-B^O. 



(3.7) 




r 
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4 Numerical Results 

The parameters corresponding to the numerical models are summarized in Table 
1. The initial flow type, the grid size, the outer radius, the Mach number at the outer 
edge Mout: the two constants of 7 and B, and the vertical thickness of the initial 
incoming flow are listed. In all these models, the inner radius is fixed at 3Rg and the 
polytropic constant is K — 4.5 x 10^^. Also note that the Mach number is defined 
excluding the rotational velocity. The Mach number is a calculated quantity in all 
models except in the Parallel- Type models in which it is a free parameter. 

The numerical code was tested in two cases in which analytical solutions exist 
relevant to the proposed study. For the case of spherical Bondi accretion (Model 1), 
the exact solution is used as the initial condition. The flow remained spherical during 
the entire calculation for each of the models tested (with varying B and K). A 
small perturbation arises at the inner edge due to the imposed boundary condition; 
however, the perturbation propagated outwards and smoothed out very quickly (in a 
timescale much shorter than the total evolution considered) . The final state is steady, 
and the relative errors of the physical parameters are always less than a few percent. 

In the second test, the evolution of a static torus (Model 2) was examined. The 
analytical solution discussed above is used as the initial condition. During the evo- 
lution, motions developed near the very surface of the torus, however, the motion in 
the interior of the torus is very small with a Mach number (excluding the rotational 
velocity) near or less than 10~^. In this case, virtually no accretion occurred. 

In the following, three distinct cases characterized by the initial conditions of the 
flow and the outer boundary conditions imposed are investigated. Specifically, they 
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are described as follows. 

Case 1: The incoming flow is supersonic and is spherically distributed with the specific 
angular momentum of gas varying with the polar angle as i{r,9) = io{l — cos 9). 
Case 2: The incoming flow is supersonic and is spherically distributed, but the specific 
angular momentum of the gas is a constant. 

Case 3: The incoming flow is supersonic and is parallel to the equatorial plane. In 
this case the thickness of the flow is restricted and the speciflc angular momentum of 
gas is a constant. 

Model 3 for Case 1 (see Table 1) is similar to the study first performed by Clarke 
et al. (1985). Qualitatively, their results that accretion occurs primarily in the polar 
direction and that a shock wave forms in the inner region near the inner edge which 
then propagates outwards is confirmed. In the subsonic regions (behind the shock) 
meridional internal motions develop and the density distribution achieves a nearly 
steady pattern. The resolution of the present calculation is higher than in Clarke et 
al. (1985) and, hence, the circulatory motion is more apparent. A typical fiow pattern 
is illustrated in Figure la on a small scale (r ^ 30Rg) and in Figure lb on a larger 
scale (r ^ 150Rg)- Here, the flow vectors correspond to a quantity related to the local 
mass fiow rate, r^pv (excluding the rotational velocity) and the contours correspond 
to the mass density. At a dimensionless time of 4830 it can be seen in Fig. la that 
the flow is primarily radial and directed toward the origin for ^ < 50 and that several 
circulatory cells are present near the mid-plane of the fiow. From Fig. lb, it is evident 
that the shock has propagated outward to r ~ 105 along the mid-plane and to r ~ 90 
in the polar direction. This difference directly reflects the presence of the centrifugal 
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forces acting only in the direction perpendicular to the polar axis. In addition, it can 
also be seen that there is a small region in the mid-plane (50 < r < 100) in which 
matter is directed away from the plane. To follow the propagation of the shock wave, 
the Mach number (vr/cs) distribution (at the mid- plane) is plotted with respect to 
the radius at four different time intervals (see Figure 2). The shock front is well 
resolved and propagates at an average dimensionless speed of 0.15, 0.075, 0.05, and 
0.025 at the four intervals indicated. It is evident that the shock decelerates as a 
function of time. We note that although the initial conditions adopted for the flow 
are not identical to that in Clarke et al (1985), the results reveal that the evolutionary 
pattern is very nearly the same. 

The primary result that an oblique travelling shock wave forms resulting in the 
development of vorticity behind the shock does not change even when the angular 
momentum distribution is constant with respect to the polar angle (Case 2). An 
example of the flow morphology is shown for Model 4 (see Fig. 3). The flow in the 
mid-plane is similar to Model 3, but the flow away from the mid-plane is significantly 
modified. In particular, there are two main differences in the flow morphology re- 
sulting from the larger amount of angular momentum in the flow. Namely, there is 
no accretion at all due to the greater contribution of the centrifugal force. That is, 
accretion in the polar direction is absent inside the funnel wall. In this wind is 

found to form near the funnel wall (r cos 9 ~ 60). However, due to the pressure of the 
incoming flow, the wind is relatively weak in this case. In fact, some of the material 
turns back to form a vortex-like flow which moves outward as a function of time. The 
development of a strong wind is favored, on the other hand, if the incoming flow is 
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artificially truncated above a given vertical height. This tendency is illustrated in 
Model 5, which is characterized by the same parameters as Model 4 except that the 
flow is truncated above H ~ Tout/ 8- It can be seen that the wind develops closer to 
the mid-plane at r cos ^ ~ 10 and extends to larger distances from the inner edge (see 
Fig. 4). 

The radial speed of the shock differs for each case reflecting the effect of the 
pressure of the incoming flow above the mid-plane of the disk. For example, the 
shock wave moves slowest in Model 5 where there is no incoming flow from above 
H «i 19. The angular momentum distribution in the incident flow can also affect the 
speed of shock propagation since the shock is found to travel faster in Model 4 than 
in Model 3. Here, the greater centrifugal support associated with the greater angular 
momentum in the flow acts in the same way to assist the enhancement of the shock 
propagation speed. 

The above simulations show that no standing shock waves are formed in the inner 
region of the disk. Even in the case of parallel incident flow characterized by a vertical 
thickness of H, standing shock waves may not form. In the following, we examine 
cases when a standing shock can be formed. Our numerical results show that the 
transition to this type of solution occurs at ~ 10, or H/vout ~ 1/3 (Model 6). 
The flow pattern for Model 6 is illustrated in Figure 5 (at time 512) where it can 
be seen that the matter immediately behind the shock front is driven out as a wind 
away from the mid-planc. No standing shock is formed inside r = 31, although it is 
possible that a standing shock exists exterior to this radius. Model 7 is corresponds 
to a boundary condition given by a parallel incident flow characterized by a vertical 
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thickness oi H — 8Rg- The flow pattern at time 4443 is shown in Figure 6. The 
sonic surface at times 4443, 2759, and 1615 is represented by the heavy sohd, dashed, 
and dotted hnes respectively. It can be seen that the shock front is approximately 
coincident with the sonic surface near r — 21 and that it is fairly steady over a long 
time interval. The shock is approximately perpendicular to the equatorial plane and 
extends in the vertical direction to about ~ 9Rg- In comparison with Model 6 matter 
in a larger region behind the shock is involved in the wind outflow. 

As the vertical extent of the incident flow is decreased further the position of the 
"standing" shock decreases. This is seen from Table 1 that the shock location is a 
sensitive function of H. In particular, the radial shock position decreases from 21Rg 
to ~ 4i?G as H is decreased by a factor of 2. The flow for a very thin incident flow 
characterized hy H = 4i?G (Model 10) is shown in Fig. 7. The shock is located in the 
very inner region of the flow at ~ 5 and rcos9 ~ 1 — 2. It can also be seen that 
the direction of the bulk flow in the wind region has also been affected by the height 
of the incident flow. Specifically, the matter is ejected closer to the "disk" surface at 
^ > 60 in comparison to Model 7 where the wind was not as restricted to he as close 
to the "disk" surface (at an angle defined by ^ > 40). 

We point out that the shock does not reach a steady state shape and position as 
the height of the incident fiow is decreased. The shock wave in Model 10 is found 
to vary in the interval 5Rg < r < SRq on a timescale of 2000, which for a compact 
object of IOMq is about 0.2 seconds. This is similar to the behavior described by Ryu 
et al (1995) for simulations characterized by incident fiows of small vertical extent. 

The effect of numerical resolution is studied in Model 11 where the number of 
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grid zones are doubled in each direction. The flow morphology is illustrated in Figure 
8 where the structure in the inner region exhibits a more well deflned vortex flow. 
However, it can be seen from Table 1 that the varibility of the shock position is 
similar to that found in Model 10, and hence the low resolution of Model 10 provides 
a reliable description of the shock variation. 

5 Discussion and Conclusion 

We have systematically investigated the structure and evolution of accretion flows 
characterized by low specific angular momentum in the adiabatic and inviscid approx- 
imation. The two dimensional hydrodynamic simulations indicate that the accretion 
behavior is rich and is primarily dependent on the distribution of angular momen- 
tum and the vertical extent of the incoming flow. For a specific angular momentum 
distribution which significantly varies with polar angle, accretion occurs in the po- 
lar direction only. Accretion in the equatorial plane is forbidden by the centrifugal 
barrier in the Newtonian approximation when the central object is located inside the 
funnel wall (Ryu et al. 1995). 

Shock waves are a common feature in these flows. For an incident flow which is 
approximately spherically distributed, the shock propagates outwards leaving behind 
a hot torus-like configuration. On the other hand, for incident fiows characterized 
by plane parallel geometry (H < IORg or H/vout < 1/3) a steady standing shock 
forms. For flows highly concentrated toward the equatorial plane {H < 5Rg) the 
shock exhibits non steady behavior. 
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In contrast to the incident flows described by a nearly spherical distribution, 
the incident flows characterized by plane parallel conditions result in the ejection of 
matter as a weak wind emanating from the inner disk region. As the flow becomes 
more highly concentrated toward the mid-plane the wind increases in strength and is 
directed away from the disk at larger polar angles. 

From this study one can conclude that the travelling shocks are transitory fea- 
tures which eventually terminate at large radii. The flow readjusts in this manner 
to the imposed initial and outer boundary conditions as it evolves to a quasi steady 
hot torus configuration. Standing shock fronts, on the other hand, are not transitory 
features in the flow. However, the signiflcance of axisymmetric standing shocks in the 
accretion process may be questioned since their applicability rests on the requirement 
that the incident flow is characterized by low specific angular momentum and a re- 
stricted spatial vertical extent. These two restrictions are usually not satisfied in most 
astrophysical applications since geometrically thin accretion usually corresponds to a 
more nearly Keplerian flow with a high speciflc angular momentum. Similarly, low 
speciflc angular momentum flows are primarily associated with geometrically thick 
flows which do not lead to standing shock structures. Only for geometrically thin 
supersonic flows where there is insufficient time to maintain a hydrostatic balance in 
the vertical direction could such structures play a role in the accretion process, e.g., 
in the case of accretion resulting from tidal disruption of a star close to a massive 
black hole. 
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Table 1 

Steady State Models 



Model Initial Flow Ne x Nr Tout Mout 7 ^/(^ H rshock 



1 


Bondi i = 


40 X 100 


149 




1.5 








2 


Torus e^io 


40 X 100 


149 





1.5 


10-5 






3 


case 1 


40 X 100 


149 


3.46 


1.5 


10-5 






4 


case 2 


40 X 100 


149 


3.46 


1.5 


10-5 






5 


case 2 


40 X 100 


149 


3.46 


1.5 


10^5 


19 




6 


case 3 


40 X 60 


31 


10 


5/3 





10 




7 


case 3 


40 X 60 


31 


10 


5/3 





8 


21 


8 


case 3 


40 X 60 


31 


10 


5/3 





6 


19 


9 


case 3 


40 X 60 


31 


10 


5/3 





5 


7-8 


10 


case 3 


40 X 60 


31 


10 


5/3 





4 


5-8 


11 


case 3 


80 X 120 


31 


10 


5/3 





4 


5-9 



21 



FIGURE CAPTIONS 

Figure 1. The flow pattern of model sequence 3 (see Table 1) at moment t = 4830. 
Panel (a) corresponds to the inner region of the flow and panel (b) to the entire region 
of the flow. In (a) the density ranges from logp = —9.12 to -7.36 (Alogp = 0.1) 
whereas in (b) the density ranges from logp = —10.52 to -7.56 (A logp = 0.1). The 
vector arrow is defined as r^pv, where v is the velocity excluding v^. The maximum 
amphtude of the vector corresponds to 4 x 10^^ g s~^. Note the development of 
vorticies in the subsonic region behind the shock and the accretion of matter near the 
polar direction. 

Figure 2. The spatial variation of the Mach number of the flow in the equatorial 
plane for model 3 for 4 different epochs. The solid, dotted, short dashed, and long 
dashed lines correspond to 472, 1412, 2978, and 4830 respectively. Note that the 
shock decelerates as it propagates outward. 

Figure 3. The flow pattern of model sequence 4 (see Table 1) at moment t — 3104. 
The density ranges from logp = —12 to -8.28 (Alogp = 0.1). The maximum vector 
corresponds to 4 x 10^^ g s~^. The development of vorticity in the inner regions is 
similar to that illustrated in Figure lb (of model 3). Here no accretion occurs as a 
result of the centrifugal barrier. An outflowing wind from the inner region is present. 

Figure 4. The flow pattern of model sequence 5 (see Table 1) at moment t — 2766. 
The density ranges from logp = —12 to -8.77 (Alogp = 0.1). The maximum vector 
corresponds to 4 x 10^^ g s~^. In this case, the wind is stronger than that in Figure 
3. In addition, the shock wave moves slower. 
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Figure 5. The flow pattern of model sequence 6 (see Table 1) at moment t — 512. 
The density ranges from logp — —12 to -7.22 (Alogp = 0.1). The vector arrow- 
corresponds to the momentum density in the flow with the largest vector equal to 
100 g cm~^ s""^. The heavy solid lines denotes a Mach number of unity. 

Figure 6. The flow pattern of model sequence 7 (see Table 1) at moment t = 4443. 
The heavy solid, dashed, and dotted lines are the sonic surface at moments t — 4443, 
2759, and 1615 respectively. The shock front which is approximately perpendicular to 
the equatorial plane at r ~ 21 is quasi-steady. Note that the vector arrow is defined 
as pv with a maximum of 70 g cm~^ s~^. The density ranges from logp = —12 to 
-7.03 (Alogp = 0.1). 

Figure 7. The flow pattern of model sequence 10 (see Table 1) at moment t = 4161. 
The heavy solid, dashed, and dotted lines are the sonic surface at moments i = 4161, 
3604, and 2041 respectively. The shock front which is approximately perpendicular 
to the equatorial plane at r ~ 5 — 8 is non-steady. Note that the vector arrow is 
deflned as the momentum density with a maximum value of 500 g cm~^ s~^. The 
density ranges from log p — —12.76 to -6.62 (Alogp = 0.1) 

Figure 8. The flow pattern of model sequence 11 (see Table 1) at moment t = 2288. 
The vector arrow is deflned as the momentum density with a maximum value of 700 
g cm~^ s~^. The density ranges from logp = —12.0 to -6.64 (Alogp = 0.1) Note that 
the vortex flow in the inner region is well resolved. 
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